Multiferroic BiFeOa-BiMnOs Nanocheckerboard From First Principles 
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We present a first principles study of an unusual heterostructure, an atomic-scale checkerboard 
of BiFeOs-BiMnOs, and compare its properties to the two bulk constituent materials, BiFeOa and 
BiMnOs. The "nanocheckerboard" is found to have a multiferroic ground state with the desired 
properties of each constituent: polar and ferrimagnetic due to BiFeOs and BiMnOs, respectively. 
The effect of B-site cation ordering on magnetic ordering in the BiFeOs-BiMnOs system is studied. 
The checkerboard geometry is seen to give rise to a a novel magnetostructural effect that is neither 
present in the bulk constituent materials, nor in the layered BiFeOa-BiMnOs superlattice. 



I. INTRODUCTION 

Artificially structured oxides present intriguing oppor- 
tunities for material design. With dramatic advances in 
epitaxial growth techniques allowing atomic-scale con- 
trol, experimental and theoretical attention has focused 
on strained-layer superlatticesi"— . Properties signifi- 
cantly different from those in the bulk have been ob- 
served, leading to the possibility of designing new mate- 
rials at the nanoscale with enhanced functionalities—"—. 
Recently, progress has been reported in the synthesis of 
artificially structured oxides with lateral "nanochecker- 
board" (or nanopillar) patterning. In particular, the 
length scale of this checkerboard ordering can be con- 
trolled by synthetic processes and stoichiometry, offering 
promise for applications such as ultrahigh-density mag- 
netic recording mediapi^"— . 

One functionality of particular current interest is mul- 
tiferroicity: the combination of ferromagnetism and fer- 
roelectricity, with coupling between the spontaneous po- 
larization and the magnetization. Room temperature 
multiferroic materials with high magnetoelectric cou- 
plings are desirable, because they can support novel func- 
tionalities in electronic devices^Si^. Magnetostructural 
and magnetoelectric couplings have been observed in a 
number of materials, including bulk^— and layered^! 
manganites, epitaxial EuTiOa^, EuSe/PbSei_a;Tea: mul- 
tilayers^^, and SrRuOa/SrTiOa oxide interfaces^!. 

Because of the distinct natures of ferroelectric and 
ferromagnetic ordering, it has proved difficult to find a 
single-phase room temperature multiferroic material with 
large polarization, large magnetization, and large mag- 
netoelectric and/or magnetoelectric coupling^^. Most 
current multiferroic devices are based on nanocompos- 
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and advances in the synthesis of artificially 
structured materials further support studies of novel 
multiferroic heterostructurcs^S. Exploring the coupling 
of ferroelectric and/or magnetic states to strain has 
shown to be exceptionally fruitful in many multiferroic 
nanocomposites^S. The challenge is to anticipate what 
new properties can arise in such heterostructures from 
combining two distinct materials, and how these proper- 
ties depend on the geometry of the combination. 

First-principles approaches are ideally suited for meet- 
ing this challenge. These methods allow searching over a 



variety of compositions, heterostructure geometries, and 
structure types to find a material with the desired prop- 
erties^. With first-principles methods, it is possible also 
to identify and characterize low-energy alternative struc- 
tures; though they are not manifest in the bulk, they 
can become physically relevant with changes in the ex- 
ternal parameters and boundary conditions produced in 
a nanocomposite. 

In this paper, we use first principles calculations to 
explore the structure and properties of a prototypi- 
cal atomic-scale checkerboard of BiPeOs and BiMnOa 
(Fig[T]), extending a shorter study of this nanocomposite 
that has been published elsewhere^. Ferroelectric anti- 
ferromagnetic (AFM) bulk BiFeOa and half-metallic fer- 
romagnetic (FM) bulk BiMnOa are good candidates for 
a nanocomposite with ferroelectric-ferromagnetic (mul- 
tiferroic) behavior. The properties of the atomic-scale 
checkerboard are found to be directly related to the prop- 
erties of the bulk constituents in their ground states and 
in low-energy alternative structures. The ground state of 
the BiFeOa-BiMnOa atomic-scale checkerboard is multi- 
ferroic, i.e. ferroelectric and ferrimagnetic, acquiring the 
desired properties from the constituents. In addition, 
we show that the BiPeOa-BiMnOa atomic scale checker- 
board displays a magnetostructural effect, namely, it 
changes its magnetic ordering with the change of its crys- 
tal structure. This effect is argued to be inherent to B- 
site cation checkerboard geometry, resulting from mag- 
netic frustration for the particular arrangement of cations 
and bonds. 

The organization of this paper follows. In See. HTJ we 
describe the first-principles method and the structural 
distortions and magnetic orderings considered. Results 
for low-energy alternative structures of bulk BiPeOa and 
BiMnOa are reported in Sees. IIIII and IIVI respectively. 
The ground state of the BiPeOa-BiMnOa atomic-scale 
checkerboard is shown to be ferroelectric and ferrimag- 
netic in Sec. |Vl A simple Heisenberg model is construc- 
tured to represent the energies of various magnetic states 
of this checkerboard computed from first principles. In 
Sec. IVIl the effect of structural distortions on the mag- 
netic ordering of the nanocheckerboard is explored, and 
we relate the properties of alternative low-energy struc- 
tures of the checkerboard to those of bulk BiPeOa and 
BiMnOa • Anisotropic epitaxial strain is shown to drive a 
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magnetic transition in the atomic-scale checkerboard in 
Sec. I VIII We study the effect of B-site cation arrange- 
ment on magnetic properties of the BiFeOa-BiMnOs sys- 
tem in Sec. I Villi The possibihty of experimentaUy re- 
ahzing a BiFeOa-BiMnOa nanocheckerboard is discussed 
in Sec. IIXI Conclusions are presented in Sec. |Xl 



alternating along z, and a rocksalt structure, with Fe and 
Mn alternating in every other unit cell ((111) superlat- 
tice in the atomic-scale limit considered here). In both 
cases, the supercell is \/2a x \pla x 2c. For consistency, 
we take the supercell for bulk BiFeOa and bulk BiMnOa 
calculations to be \f2a x \/2a x 2c, except for the R3c 
structure, where we use a \f2a x \f2a x \/2a supercell. 



II. METHOD 

First-principles calculations are performed using den- 
sity functional theory within the local spin-density ap- 
proximation with Hubbard U (LSDA-I-U) method as 
implemented in the Vienna ah initio simulation pack- 
age VASP-4.6.34^2i^. Projector-augmented wave po- 
tentials (PAW)^i^ are used with 15 valence electrons 
for Bi (5di°6s26p3), 14 for Fe (Sp^Sifis^), 13 for Mn 
(Sp^Sd^As"^), and 6 for O {2s^2p^). The robustness of 
the results is tested with two different implementations 
of the rotationally invariant LSDA+U version. The first 
is due to Liechtenstein^ with effective on-site Coulomb 
interaction Upe = C^j\/n = SieV and effective on-site ex- 
change interaction Jpe = Jmu ~ leV. The second is 
due to Dudare\^, with Ul{l = 5.2eV, VjJ = AeV, 
where If^^-^ = U — J. Both implementations treat lo- 
calized d electron states in Fe and Mn. It has previ- 
ously been shown that these U and J values give good 
agreement with experiment in bulk BiFeOa^. The value 



U 



Mn 



5.2eV has previously been used for bulk BiMnOa 



ground state calculation: 
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Figure 2: Magnetic orderings considered for bulk BiFeOa 
and bulk BiMnOs: Symbols for each type of ordering are 
introduced next to each label. 



Fe|-Mnf-7| 



Fe|..Mn|.^/^ 



G-AFM O 



Fe|--Mn|l--7| 
Fe|_.Mn|.^ 



C-FIM □ 



Fe:| 



Sot 



Fe|..MnJ^/^ 



FM A 



FeJ-Mnf-Ti Fe|--Mn|--7| FeJ--MnJ--7| 

nw\ 

FeJ..Mn|.^ FeJ_.Mn|.^ FeJ_.Mn|.^ 



FeAFMMnFM ♦ FeFMMnAFM ▼ FeAFMMnAFMO 



Figure 3: Magnetic orderings considered for the BiFeOs- 
BiMnOs atomic-scale checkerboard compatible with the 
y/2a x \/2a x 2c supercell. Symbols for each type of ordering 
are introduced next to each label. 



Figure 1: (a) BiFeOa-BiMnOa atomic-scale checkerboard, 
(b) Top view of the atomic-scale checkerboard, (c) Perovskite 
cell. Dashed lines show an oxygen octahedron surrounding the 
B-site (Fe, or Mn) cation. 



The BiFeOa-BiMnOa atomic-scale checkerboard is 
shown in Fig. [T] BiFeOa and BiMnOa alternate at the 
atomic level, forming a checkerboard pattern in the xy 
plane and pillars of the same composition along z. The 
supercell is \/2a x y/2a x 2c, containing two Fe and two 
Mn. In the limit of the atomic-scale pillars considered 
here, the checkerboard structure is the same as that of a 
(llO)-oriented superlattice. 

We consider two additional types of B-site cation- 
ordered BiFeOa-BiMnOa systems: a (OOl)-oriented lay- 
ered superlattice, with single unit-cell Fe and Mn layers 



Several types of magnetic orderings are studied here: 
the G-type (rocksalt), C-type, A-type AFM, and FM or- 
dering of the local magnetic moments in bulk BiFeOa, or 
bulk BiMnOa (see Fig. [2]). All orderings considered are 
collinear; this is supported by recent neutron scattering 
measurements on BiFeOa doped with Mn^, that indicate 
collinear AFM ordering. 

Magnetic orderings of the BiFeOa-BiMnOa atomic- 
scale checkerboard are shown in Fig.|3l where we consider 
six collinear orderings of Fe and Mn spins. Similarly, six 
collinear orderings of the magnetic Fe and Mn spins are 
explored in the BiFeOa-BiMnOa (OOl)-oriented superlat- 
tice and the rocksalt structure. For the (OOl)-oriented 
superlattice, these orderings are described by the nota- 
tion FcFM (FeAFM), or MnFM (MnAFM), referring to 
the FM (AFM) ordering for the Fe (Mn) moments in 
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the relevant layer, respectively, with the remaining am- 
biguities resolved as follows: FeAFMMnAFM magnetic 
order has AFM ordered Fe and Mn layers with FM order 
along the mixed Fe-Mn chains in the z direction, while 
G-AFM designates the case with AFM order along the 
mixed chains; similarly, FeFMMnFM has FM ordered Fe 
and Mn layers with AFM order, while FM designates 
the case with FM order along the mixed chains. For 
the rocksalt structure, we consider FM and G-AFM or- 
dering, FeAFMMnFM ordering, referring to FM ordered 
Mn sublattice and AFM ordered Fe sublattice; similarly 
we consider FeFMMnAFM ordering with FM ordered Fe 
and AFM ordered Mn sublattices, respectively. Finally, 
the FMFM ordering has AFM ordered Mn and Fe sub- 
lattices, which are coupled FM in each Fe-Mn z layer, 
while the AFMAFM ordering has AFM ordered Mn and 
Fe sublattices coupled AFM in each Fe-Mn z layer. 
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Figure 4: Distortions of the ideal cubic perovskite cell con- 
sidered in this work^: (a) Polar distortion with relative shift 
of Bi and Fe/Mn cations with respect to O anions along the 
2-axis (rj(z)), (b) + (in-phase) rotations of the oxygen octa- 
hedra (dashed lines)about the z-axis [M'^(zy), (c) — (out-of- 
phase) rotations of the oxygen octahedra (dashed lines) about 
the z-axis (_RJ(2)). 



Structures generated by three modes of the cubic per- 
ovskite structure are considered (see Fig. U])^: (i) the 
zone center polar F^ mode, (ii) the oxygen octa- 
hedron rotations (all rotations about a given axis are in 
phase), and (iii) i?^ rotations (sense of rotations alter- 
nates along the rotation axis) . Space groups correspond- 
ing to the structural distortions considered are presented 
in Table HI and we use the notation c-, 1-, or r- to refer 
to the structural distortion of the B-site cation-ordered 
checkerboard, layered superlattice, or rocksalt structure, 
respectively. To search for the ground state for a given 
magnetic ordering and structure type, we perform struc- 
tural relaxation with the conjugate gradient algorithm. 
Both the cell shape and the cell volume are relaxed; more 
specifically, the ions are relaxed towards equilibrium po- 
sitions until the Hellmann-Feynman forces are less than 
10~'^eF/A. An energy cutoff SSOeF for the plane wave 
basis set is used. Convergence in the energy is reached 
with precision 10~ ''eF. A Monkhorst-Pack k-point grid^ 
is generated with density 4 X 4 X 4 for the (%/2 X ^/2 X ^2) 



supereell, and 4 x 4 x 2 for the (v^ x \/2 x 2) supercell. 
For magnetic energy calculations (Sees. IVB] and IVlIIl) . 
we use the energy cutoff SOOeT^, and the Monkhorst-Pack 
k-point grid with density 6x6x4. Gaussian broad- 
ening of the partial occupancies for each wavefunction 
is 0.05eF. A tetrahedron method with Blochl correc- 
tions^ is used for the density of states (DOS) calcula- 
tions, with the Monkhorst-Pack k-point grid 4 x 4 x 4 for 
the (V2 X V2 X \/2) and 8 X 8 X 4 for the (\/2 x V2 x 2) 
supercell. 

The rotational distortion can be quantified using the 
oxygen octahedron rotational angle defined specifically 
for each oxygen in the octahedron as 



COS0 



(1) 



where u is the shortest vector from the rotation axis to 
the reference position of the oxygen, and v is the shortest 
vector from the rotation axis to the position of the oxygen 
in the distorted structure. The rotation axis is [001] and 
[010] for the M;^(z) and R^iy) distortions respectively, 
and the threefold axis (body diagonal of the cube or dis- 
torted cube) for the i?J([lll]) distortion. Due to defor- 
mation of the oxygen octahedron in the BiFeOs-BiMnOs 
checkerboard structures, these angles may be different for 
different oxygens in the same octahedron. We report an 
average value if the range is small; otherwise the lower 
and upper limits of the range are presented. 

The polar distortions of the various structures can be 
quantified by estimating the polarization based on a lin- 
earized expression with nominal charges: 



'J' 



(2) 



where P is the polarization, Auj is the displacement of 
the jth ion with respect to its ideal perovskite position, qj 
is the nominal charge of the jth ion (qBi = +3, qpe = +3, 
QAin = +3, qo = ^2), and is the unit cell volume. 

For selected structures the true value for the sponta- 
neous polarization is computed using the Berry phase 
method^i^ as implemented in VASP-4.6.34. In this for- 
malism, the polarization is only well-defined mod ei?/f2, 
where R is any lattice vector and is the primitive- 
cell volume; thus possible values of the polarization are 
points on the lattice defined by Pq + eR/il, where Pq is 
the value directly obtained from the Berry phase calcu- 
lation. Choosing the lattice point (or branch) that corre- 
sponds to the measured switching polarization (e.g. in an 
electrical hysteresis loop) is done by computing the po- 
larization of states closely spaced along an adiabatic path 
connecting the structure of interest to a high-symmetry 
reference structure. These laborious calculations can be 
avoided by an approach based on the reformulation of 
the polarization in terms of Wannier function centers^; 
the switching polarization is obtained from the differ- 
ence between the two symmetry-related variants by as- 
sociating the Wannier centers with the same atoms in 
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Table I: Resulting space groups for considered structural distortions (see Fig. |4)|. Pm3m is the ideal perovskite structure and 
P4/mmm is the uniformly strained tetragonal unit cell. 



Modes r^{z) M+iz) M+(z),rj(z) R+{y) i?+(y),r,-(y) i?+([lll]) 7?+([lll]),rr ([111]) 
Abbrevation (z) (z) Mr{z) Rt{y) Rr(y) Rt{d) Rr{d) 

Space Group P4mm P4/mbm P4bm I4/mcm I4cm R3c R3c 



both structure a ' . Due to incompatibility between the 
WannierQO and the VASP codes, we cannot use this lat- 
ter approach here; we make the necessary branch choices 
based on computations along adiabatic paths combined 
with the nominal-charge polarization estimate. 



Fe d states. The polar character of BiFcOa arises from 
the polar rj([lll]) mode, and the spontaneous polariza- 
tion using the Berry phase method is P^-^c _ gO/iC/cm^ 
along the [111] direction. 



III. BIFEO3 STRUCTURES 
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Figure 5: Density of states (DOS) of the R3c G-AFM ground 
state (GS) of BiFeOa. Spin up states are plotted by solid 
line and spin down states by a dashed line (inverted view). 
The zero is set to the valence band maximum. DOS shows an 
insulating band gap of 2 eV. 



In agreement with previous first-principles calculations 
and experiment^i^i^i^, we find that the ground state 
structure of BiFeOa has rhombohedral R3c symmetry, 
which is a combination of the rotational i?^([lll]) mode 
(counter-rotations of the oxygen octahedra about the 
[111] axis) and a polar r4 ([111]) modes, with Bi, Fe, 
and O displaced relative to one another along [111] and 
further distortion of the oxygen octahedra by displace- 
ment of the O displaced perpendicular to [lll]4i The 
oxygen octahedra rotation angle is large, about 14°, and 
is comparable to rotations in other perovskites^. The 
R3c ground state has G-AFM (rocksalt) ordering (see 
Fig. [5]), and Fe local magnetic moment of A^b^ The 
density of states (DOS) is plotted in Fig. [5J it has a 
2eV band gap that separates occupied and unoccupied 
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Figure 6: Structural energetics of bulk BiFeOg. Energy differ- 
ence per Fe for different magnetic orderings (see Fig. [2]) and 
structural distortions (see Fig. |4] and Table [l| relative to the 
FM Pm3m structure. 



Next alternative structures of BiFeOa are studied, and 
we consider those generated by freezing in linear combi- 
nations of the rotational ^ , and polar F^ modes, 
and four magnetic orderings (see Fig. [2]) . Their energies 
are plotted in Fig. [H] (see also Table HIl) . By symmetry, 
the FM ordering has the same energy for the (z) and 
F^ (y) structures; this is also true for the G-AFM order- 
ing. In contrast, the C-AFM ordering has different en- 
ergy for the T^{z) and r4 (j/) structures; this is also true 
for the A-AFM ordering, as the y and z directions for 
these spin arrangements are not symmetry-related (see 
Figs. [2 and [H). 

For all structural distortions considered, the favored 
magnetic ordering is G-AFM (open circle). This is con- 
sistent with the Goodenough-Kanamori rules: either a 
strong (T bond is formed between Fe eg and the neigh- 
bouring O p orbitals in an ideal 180° Fe-O-Fe bond (ideal 
perovskite structure), or a weak tt bond is formed be- 
tween Fe t2g and O p orbitals when the bond is bent 
towards 90° (as the structure is distorted); in both cases, 
the AFM superexchange is favored^"—. 
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Table II: GS and low-energy alternative structures of bulk 
BiFeOs. The energy difference AE is given with respect to the 
FM Pm3m structure, as in Fig.[6l Polarization P is estimated 
from the nominal charges (Eq. 0). Also included are the 
band gap A and the a and c lattice constants for the V2a x 
y/2a X 2c supercell of P4mm, and the a lattice constant and 
the angle a for the \/2a x \/2a x y/2a supercell of R3c. 



Space group 


P4mm 




R3c 


Modes 






i?+([iii]), r4-([iii]) 


Mag. order 


G-AFM 


C-AFM 


G-AFM 


A£;[eV/Fe] 


-1.09 


-1.10 


-1.25 


A [eV] 


1.75 


2.23 


1.99 


P[fj.C/cm'^] 


113.6 


116.2 


62.1 


a/ c [A] 


3.68/4.64 3.67/4.68 


5.52, Q = 59.8° 



IV. BIMNO3 STRUCTURES 

Previous first principles calculations show that the 
ground state structure of bulk BiMnOs is monoclinic 
centrosymmetric C2/c with zero spontaneous polariza- 
tion^ and FM ordering—"—. Bismuth (Bi) cations are 
ofF-ccntcr due to stcreochcmically active Bi lone pairs, 
and the Jahn- Teller activity of Mn'^^ further distorts the 
structure^. Optimizing the atomic positions and lattice 
constants, we performed a first-principles calculation for 
this structure to find an energy gain of 1.26eV /Mn rel- 
ative to the the ideal cubic perovskite structure with G- 
AFM ordering and qq = 3.83A; the latter is used as our 
reference state throughout this section. 

We study low-energy alternative structures of BiMnOa 
compatible with a V2 x \/2 x 2 supercell. Results are 
presented in Fig. [H The lowest energy structure has R3c 
symmetry, the same structure type as the ground state of 
BiFcOa. It is FM, with magnetic moment S.QfiB per Mn. 
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The most favorable low-energy alternative structures 
and the ground state of BiFeOa are presented in Ta- 
bleini The low-energy (z) structure, with P4mm sym- 
metry, is the supertetragonal structure with c/a 1.3, 
previously discussed elsewhere^^^. It has been shown 
recently that this phase can be stabilized in BiFeOs thin 
films^. The polarization computed by the Berry phase 
method^ is p^^mm ^ 150/iC/cm^. The nominal-charge 
estimates are therefore smaller than the true values for 
both the P4mm and R3c structures, but the relative val- 
ues are well reproduced. The polar distortion is dominant 
in all structures considered; the rotation-only structures 
are higher in energy, and the presence of a polar distor- 
tion tends to stabilize the rotational instabilities. For ex- 
ample, the -I-A/3'' and T^+R^ distortions relax back 
to with zero oxygen octahedron rotation angle. The ^ 
only exception is the R3c structure, in which the rota- 
tional and polar distortions coexist. c 
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Figure 7: Density of states (DOS) of the R3c FM alternative 
structure of BiMnOa. Spin up states are plotted by solid line 
and spin down states by a dashed line (inverted view). Fermi 
energy is shown by the vertical dashed line crossing zero. DOS 
shows a half-metallic character. 



This structure lies only 43mey/Mn above the BiMnOa 
monoclinic ground state. The computed DOS is shown 
in Fig. [7j the system is half-metallic, with a gap of 3.25 
eV in the spin down channel. As an aside, we note that 
it might be useful to stabilize BiMnOa as a half-metal 
in this low energy structure for possible applications in 
spintronics^. 
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Figure 8: Structural energetics of bulk BiMnOs. Energy 
difference per perovskite cell (Mn) for different magnetic or- 
derings (see Fig. [2]) and for structural distortions (see Fig. |4] 
and Table H]) generated by the specified modes. 



For all structural distortions considered, the favored 
magnetic ordering is FM, consistent with previous anal- 
ysis that showed that BiMnOa favors FM structures 
with a half-metallic characte r^^-^'^ . The ferromag- 
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netism in BiMnOs can be explained by a combina- 
tion of Goodenough-Kanamori rules and orbital order- 
ing^i^i^i^. Structural distortions (cither oxygen oc- 
tahedron rotations or polar distortion) widen the spin- 
down gap (see Table llll| : a similar trend is observed for 
the band gap in BiFeOa (see Table A small band gap 
opens with a monoclinic distortion in the FM BiMnOa 
ground stat o^"'^^ . 

The Jahn- Teller active Mn'^+ configuration tends to 
favor elongation of the oxygen octahedron. In contrast 
to BiFeOs, in which the polar instability strongly dom- 
inates, the rotational and polar instabilities in BiMnOa 
are comparable in magnitude, as can be seen by compar- 
ing the energies of the F^ (z), M^{z) and R^{y) states. 
The latter two states have a small residual polar instabil- 
ity. A polar distortion along a Cartesian axis lowers the 
energy of the G-AFM state so that the energy difference 
between this state and the FM ground state is greatly re- 
duced; this does not occur if the polar distortion is along 
[111] as in the R3c phase. The octahedral rotation an- 
gles in the low-energy BiMnOa structures are all similar 
in magnitude, varying between 11 — 14°, with an angle 
of 13° for the FM R3c structure. The value of the octa- 
hedral rotation angle in the G-AFM R3c structure, 14°, 
is the same as in G-AFM R3c BiFeOa. 



the net magnetization cancels due to the AFM pillar or- 
dering, the contribution from Mn moments adds, leading 
to a net magnetization of 3.8/xb per Fe-Mn pair. 
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Figure 9: Density of states (DOS) of the c-R3c FeAFMM- 
nPM ground state (GS) of BiFeOa-BiMnOa nanochecker- 
board. Spin up states are shown by a solid line and spin 
down states by a dashed line. The zero is set to the valence 
band maximum. The band gap in the spin-up channel is 0.9 
eV. 



V. BIFEO3-BIMNO3 NANOCHECKERBOARD 
GROUND STATE 

A. Crystal structure, magnetization and 
polarization 

In the search for the ground state of the atomic- 
scale checkerboard cation ordering, we considered the 
six coUinear magnetic states of Fig. [3] and four different 
structures: the tetragonal P4/mmm structure and three 
additional structures, obtained by freezing in a F4 (z) 
mode, a combination of R4^(y) and F4 (y), and a com- 
bination of R4 (111) and r4 (111). We designate these 
latter three structures by the space group they would 
have if all B sites were occupied by the same cation, with 
the prefix c- to remind us that the actual symmetry is 
lower due to the checkerboard ordering: c-P4mm, c-I4cm, 
and c-R3c. The GS of the BiFcOs-BiMnOs nanochcckcr- 
board is found to be c-R3c, as could be expected based 
on the R3c GS of bulk BiFeOa, and on our results for 
bulk BiMnOs. The magnetic ordering in the c-R3c GS is 
FeAFMMnFM. Fe magnetic moments are ordered AFM 
along the Fe pillars and Mn magnetic moments are or- 
dered FM along the Mn pillars, as expected from the G- 
AFM and FM ground states of BiFeOa and BiMnOa, re- 
spectively (see Sees, mil and lIVp . AFM and FM xy layers 
alternate along z as is sketched in Fig. [3] The computed 
Fe and Mn local magnetic moments are 4.1^b and 3.8/iB, 
respectively; these are the same values as those reported 
here in the parent compounds BiFeOa and BiMnOa. Al- 
though the contribution from Fc magnetic moments to 



The DOS of the c-R3c FeAFMMnFM GS is shown in 
Fig. ini The general features are very similar to those 
found in BiFeOa (Fig- El) and BiMnOa (Fig.©, the main 
difference being that the spin-up Mn states at the Fermi 
level in BiMnOa have split to open a gap, with the occu- 
pied states at the top of the valence band narrowing the 
gap to 0.9eV. 

Direct calculation of the spontaneous polarization us- 
ing the Berry phase method for the c-R3c GS yields a 
valueof P= (19.6, -1.1, 30.5) fiC/cm^. This is well de- 
fined only up to the polarization lattice eR/fl^^^, which 
in this case is (13.2, 13.1, 0.1)ni-h(-13.2, 13.0, 0.1)n2 + 
(0.0, 0.2, 26.6)7ia isC/cw?, where n is a vector of in- 
tegers. To determine the branch that corresponds to 
the switching polarization, we compute the polarization 
along a structural deformation path that linearly con- 
nects the c-R3c GS to the ideal cubic perovskitc struc- 
ture. As shown in Fig. 1101 the computation is performed 
for structures down to 75% of the full distortion (at which 
point the structures become metallic) and then linearly 
extrapolated to 0% using the expression 

APioo%-o% = 4xAPioo%-75% = (33.8, 33.5, 39.8)MC/cm2. 

(3) 

The magnitude of this estimate, 62.0/iC/cTO^, suggests 
the branch choice P(,p = (32.9, 38.0, 30.7) ^C/cm? with 
magnitude \Pbp\ = 58.9^C/cm^. However, we would ex- 
pect it to be considerably larger, based on comparison 
between the Berry phase |P| of bulk BiFcOa and the 
\P\ computed using nominal charges (see Sec. IIII[) : fol- 
lowing this intuition we would make the branch choice 
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Table III: Low energy alternative structures of FM bulk BiMnOs. The energy difference A_E is calculated with respect to 
the G-AFM Pm3m structure (as in Fig. Listed are values of the spin-down band gap Ahm in the half-metallic structures 
or metallic (m) character, the oxygen octahedron rotational angle Q (see Sec. |Il]| , and the a and c lattice constants of the 
-/2 X \/2 X 2 supercell. 



Space group P4/mmmm P4bm I4cm R3c 

Modes M+{z), V^{z) R+{y), F^jy) fi|([lll]), F^ ([111]) 

AE [eV/Mn] -0.57 -1.02 -1.03 -1.22 

Ahm [eV] m 0.73 2.74 3.25 

e[°] - 11.4 12.0 13.3 

a/ c [A] 3.83/ 3.86 3.81/ 4.01 3.81/ 3.83 5.51, a = 60° 
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Figure 10: Spontaneous polarization P = [Px, Py, Pz) as a 
function of the structural distortion of the c-R3c FeAFMM- 
nFM BiFeOs-BiMnOa nanocheckerboard. 100% distortion 
corresponds to the c-R3c ground state structure, and 0% to 
the ideal perovskite structure. 
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Figure 11: Sketch showing the displacements of oxygen 
atoms in the i?^([lll]) mode, which contributes to the c-R3c 
ground state (GS) of the BiFeOs-BiMnOs nanocheckerboard. 
Two inequivalent iron (Fei and Fe2) and manganese (Mui 
and Mn2) atoms and twelve oxygens (1 — 12) comprise the 
72 X 72 X 2 unit cell. The corners of each cube are occupied 
by Bi (not shown). 



Pbp = (46.1, 51.3, 57.4) fiC/cm^ with magnitude \Pbp\ = 
89 .7 iiC / cm^ . This remaining ambiguity highlights the 
challenge of picking the right branch when the polar- 
ization is much larger than the quantum; in either case 
it is clear that the polarization of the checkerboard is 
comparable to the largest values found in ferroelectrics. 
Thus, we find that the c-R3c GS of the BiFeOs-BiMnOs 
nanocheckerboard is multiferroic: ferroelectric, with po- 
larization comparable to the polarization of bulk BiFeOs, 
and ferrimagnetic, with magnetization contributed by 
ferromagnetic ordering in the BiMnOa component. 



B. Magnetic Coupling Constants 

To gain insight into the magnetic properties of the 
nanocheckerboard, we model the magnetic ordering en- 
ergies using a nearest-neighbor (nn) Heisenberg model. 



The nn magnetic couplings arise from superexchange 
through the oxygens that lie on the bonds between the 
B site cations, with the strength of the superexchange 
being quite sensitive to the geometry of the B-O-B bond. 
If the structure were ideal cubic perovskite, there would 
be three independent couplings, Jpe, ^Mn and Ji„t, cor- 
responding to 180° Fe-O-Fe, Mn-O-Mn and Fe-O-Mn 
bonds, respectively. The analysis of the couplings in the 
c-R3c structure is based on the geometry of the B-O-B 
bonds as given in Table IIVI the labeling of the bonds 
and the changes in the bonds due to the Rj([lll]) mode 
are shown in Fig. [TTJ The two Fe-O-Fe bonds are almost 
identical in bond angle and bond length; this is also the 
case for the two Mn-O-Mn bonds. This suggests that a 
single value of Jpe and Jmu can be used for the Fc-O-Fe 
and Mn-O-Mn interactions, respectively. On the other 
hand, the mixed Fe-O-Mn bonds vary in both B-0 bond 
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Table IV: B-site-cation-oxygen-B-site-cation (B-O-B) bonds in the c-R3c ground state (GS) of the BiFeOa-BiMnOs 
nanocheckerboard. B-O and 0-B bond lengths and the B-O-B bond angle are given. Atoms are numbered as in Fig. 1111 
The subscript indicates the cartesian direction along which the bond lies. 
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length, from 1.91 — 2.1lA, and B-O-B bond angle, from 
143.3 — 166.8°. This suggests the use of three different 
coupling constants J^^., jf^^ or J^^^ for the Fe-O-Mn in- 
teractions based on the typical values of the bond angles, 
approximately 166°, 154° and 144°, respectively. Note 
that the angles of the Fe-O-Mn bonds in the Fci-Mui 
layer in Fig. [Tl]are about 166° and 154°, while the an- 
gles in the Fe2-Mn2 layer are about 152° and 144°. 

The values of these five exchange couplings were deter- 
mined from first-principles results for the total energies of 
various magnetic orderings for the c-R3c GS structure of 
the nanocheckerboard, given in Table |Vl The structure 
is fixed to that obtained for the FeAFMMnFM ordering 
Fe| Fe^ Mn| Mn^. The ordering Fef Fe| Mn| Mn|, also 
described as FeAFMMnFM, is a distinct state with a dif- 
ferent (higher) energy. Similarly, for the FeFMMnAFM 
ordering, there are two distinct magnetic states: Fe| Fe| 
Mnf Mn| and Fej Fe| Mn| Mn^, with different energies 
as given in Table IVl 

We express the Heisenberg magnetic energy of each 
magnetic state, 

^^ = ^^o + i^J.,5.5„ (4) 

where Si and 5*^ are the spins Spe ~ fj^A/n = | with 

coupling constants Jy = Jpe-, Jmu, JSit^ ^Ziv ^^d 
Eq is a constant. We extract values of the coupling con- 
stants by fitting the Heisenberg model energy to the first- 
principles energies by the least-squares method, obtain- 
ing 

Eq = 218meF, Jpe = 7.1meV, = — 3.2meF, 

JS^t = ~3.0meV, jf„, = 4.3meF, Jl,^7.3meV. 

(5) 



The quality of the fit can be assessed by comparing the 
first-principles energy to the fitted values in the fifth col- 
umn of the table. 

The AFM character of Jpe and the FM character of 
Jmu correspond to that of bulk G-AFM BiFeOs and bulk 
FM BiMnOa, respectively. Their values are comparable 
to those obtained from the observed bulk transition tem- 
peratures within mean field theory assuming a single J: 
JpeMk ~ G.SmeV and JMnMik ~ -l.5meV, respec- 
tively^"—. The correspondence is not exact because of 
the difference in the bond geometry between bulk BiFeOa 
and bulk BiMnOa and the nanocheckerboard. 
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Figure 12: Angular dependence of the Fe-Mn magnetic ex- 
change coupling constant Jint in the c-R3c GS structure of the 
BiFeOa-BiMnOs nanocheckerboard (blue circles). The black 
line is a linear fit. Ji„t changes sign (AFM Ji„t > to FM 
Jint < 0) at the bond angle of 160°. 



In the Fei-Mui layer, the average Fe-Mn interaction 
J,nt is very weakly AFM {J^^^ + J^t)/'^ ^ 0- while in 
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Table V: Calculated magnetic energies /S.E per four-perovskite unit cell (u.c.) in the c-R3c GS structure of BiFeOa-BiMnOs 
nanocheckerboard. The notation for magnetic ordering is that of Fig. [S] The symbols x, y, a, b, c appearing in the magnetic 

energy are defined as follows: X = JpeSFeSpe, y = JMnSMnSun, a = J"„tSFeSM„, b = jf^^SpsSMn, C = Ji^tSpeSMn- 
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the Fe2-Mn2 layer it is strongly AFM ( 



-Jlt)/2>0- 



This corresponds to the preferred FeAFMMnFM order- 
ing Fe| Foj Mnj Mn^ and explains the close competition 
with G-AFM ordering, in which both layers are AFM or- 
dered (see Fig. [3] and Table El) . The exchange coupling 
between Fe and d'* Mn takes place via superexchange 
through the bridging O. For angles close to 180°, strong 
a bonding favors FM ordering. However, as the Fe-O-Mn 
angles deviate from 180° through the oxygen octahedron 
rotational distortion, the admixture of tt bonding leads 
to an increasingly AFM character of the coupling^"—. 
This behavior can be seen in the dependence of the fitted 
values for Jint on the Fe-O-Mn angle, plotted in Fig. [121 

Within this nearest-neighbor Heisenberg model, we ex- 
plored a wider range of possible magnetic orderings for 
the R3c structure nanocheckerboard, in particular, or- 
derings with lower translational symmetry than those in- 
cluded in the first-principles investigation. The supercells 
considered included 2x2x2 (p = 8 perovskite cells), 
2x2x4 and 4 x 2 x 2 (p = 16 perovskite unit cells). The 
Heisenberg model energies were computed for all 2^ spin 
configurations in each supercell. 

The lowest energy ordering found in this larger set of 
configurations is still the FeAFMMnFM ordering, with 
FM alignment of the Mn and AFM antialignment of the 
Fe along the Mn and Fe pillars, respectively, and alter- 
nating FM and AFM xy layers as in Fig. [31 The lowest- 
energy alternative magnetic state is a state in which one 
Mn per supercell in the FM xy layer flips, at an energy 
cost of 6.3 meV/supercell. The net magnetization for 
the resulting state decreases from Mqs = 3.8/iB per one 
Fe-Mn pair to ^Mqs- 

Within a mean field approximation with four effec- 
tive fields, two for the two Fe atoms and two for the 
two Mn atoms in the unit cell of the BiFeOa-BiMnOs 
nanocheckerboard, the magnetic transition temperature 
of the BiFeOa-BiMnOa nanocheckerboard is = i06K. 
This temperature is intermediate between the Neel tem- 
perature, = 643 K, of bulk BiFeOa and the Curie 



temperature, T^'^p = 105 K, of bulk BiMnOa^"^. 



VI. ALTERNATIVE STRUCTURES OF THE 
BIFEO3-BIMNO3 NANOCHECKERBOARD 
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Figure 13: Structural energetics of BiFeOs-BiMnOs 
nanocheckerboard. Energy difference per perovskite cell (per 
B-cation) for different magnetic orderings (see Fig. (3]) and 
for structural distortions (see Fig. [H and Table [l|. Inset: 
zoomed view of the magnetic energies for the c-R3c struc- 
ture. FeAFMMnFM ordering (filled diamond) competes with 
G-AFM ordering (open circle) in the nanocheckerboard alter- 
native structures. 



The energies for various magnetic orderings and struc- 
tural distortions of the nanocheckerboard are shown in 
in Fig. 1131 The structural parameters for each struc- 
ture type are relaxed for each magnetic ordering. The 
most energetically favorable alternative structures, like 
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Table VI: Low energy alternative and ground state (GS) structure of BiFeOa-BiMnOs nanocheckerboard. Energy difference 
A.E[eV/B — cation] is calculated for different magnetic orderings with respect to FeAFMMnAFM P4/mmm structure [as in 
Fig. [13]. Insulating DOS band gap A (or metallic m character), polarization P estimated from the nominal charges [Eq. ([2])], 
the in- plane a and out-of-plane c lattice constants [see perovskite cell in Fig. [1], and oxygen-octahedron rotation angle O are 
given for the lowest-energy magnetic ordering corresponding to each structural distortion. 



Space group 


P4/ mmm 


c-P4mm 


c-I4cm 


c-R3c 
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- 
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Mag. order 
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1.01 
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-0.913 


A£; (C-FIM) 


-0.042 


-0.654 


-0.774 


-0.884 


AE (FeFMMnAFM) 


-0.054 


-0.625 


-0.717 


-0.865 


A£ (FeAFMMnAFM) 


0.000 


-0.592 


-0.721 


-0.856 


AE (FM) 


-0.077 


-0.567 


-0.661 


-0.817 



the ground state FeAFMMnFM c-R3c structure, are po- 
lar and include oxygen octahedra rotations. 

The polar distortion in the alternative structures of the 
nanocheckerboard is quantified by the value of the polar- 
ization computed using nominal charges (Eq. that 
can be directly compared with nominal-charge polariza- 
tions in the structures of BiFeOs (cf. Tables HIl and IVTl) . 
As in BiFeOa, there is a low- lying supertetragonal P4mni 
phase, with c/a 1.3 and very large spontaneous polar- 
ization. For the various structures considered, the po- 
larization tends to decrease as rotational distortion is 
introduced, with the smallest value found in the c-R3c 
structure. 

In the P4/mmm structure, the nanocheckerboard is 
metallic, while a band gap opens with cither polar or 
rotational distortion. This behavior is similar to that 
of BiMnOa and BiFeOa, which are metalhc in the FM 
Pm3m, or P4/mmm structures with a band gap opened 
and/or widened by distortion (in FM BiMnOa, only in 
the spin-down channel). 

As can be seen in Fig.[13l the difference in energies be- 
tween different structure types is generally much larger 
than the difference in magnetic energies for a given struc- 
ture type. The interesting feature of this figure is that the 
favored magnetic ordering is different for different struc- 
ture types, switching between ferrimagnetic FeAFMM- 
nFM and antiferromagnetic G-AFM. This is in contrast 
to the case of bulk BiFeOs (see Fig. H), or bulk BiMnOa 
(see Fig. [5]) , in which the favored magnetic ordering does 
not change for different structure types. 



VII. MAGNETIC AND STRUCTURAL 
TRANSITIONS DRIVEN BY ANISOTROPIC 
STRAIN 

The sensitivity of the magnetic exchange couplings to 
the structure should produce changes in the magnetic 
ordering energies for perturbations that couple to the 
crystal structure, such as electric field, pressure and epi- 
taxial strain. It is even possible that a structural per- 
turbation could drive the system through a magnetic 
transition into an alternative low-energy ordering. Fur- 
thermore, the fact that in the nanocheckerboard the fa- 
vored magnetic ordering is different for different structure 
types, discussed in the previous section, suggests that 
the magnetic ordering of the system could in principle 
be changed by a perturbation that changes the structure 
type, producing a novel magnetic-coupling response at 
the magnetic-structural phase boundar y^^i'^^ . For exam- 
ple, it might be possible to drive the nanocheckerboard 
from its ferrimagnetic FeAFMMnFM c-R3c GS with a 
nonzero magnetization to a G-AFM c-I4cm state with 
zero magnetization. 

We have explored this possibility for two types of epi- 
taxial strain. First, we investigated the c-R3c phase with 
an isotropic epitaxial strain, corresponding to an (110) 
matching plane. Thus, the second and third lattice vec- 
tors of the y/2a x \/2a x 2c supercell, along [-110] and 
along [001], are constrained to be perpendicular with uni- 
form scaling of the lattice constants a = {1 -\- s)ao and 
c = {1 -\- s)co, where ao = 5.50A and cg = 3.93A are the 
unstrained lattice constants of FeAFMMnFM c-R3c GS. 
In this case, there is no magnetic transition: the system 
remains FeAFMMnFM from s = 0% up to strain of 10%. 

Second, we considered an anisotropic epitaxial strain. 
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Figure 14: Epitaxial-strain-driven magnetic transition in 
BiFeOs-BiMnOa nanociieckerboard. Total energies of the 
FeAFMMnFM (solid) and G-AFM (dashed line) magnetic or- 
derings in the c-R3c structure type (blue) and c-I4cm struc- 
ture type (red), as a function of anisotropic in- plane tensile 
strain (details in the text). 



corresponding to a (001) matching plane, such that 
the lattice constant along [110] is fixed to \f2 x oq = 
5.52A, while the lattice constant a along [-110] is elon- 
gated, with strain defined as . Og = 3.9A is cho- 
sen as it is the lattice constant of an ideal perovskite 
cell with volume which is the average of that of bulk 
R3c BiFeOs iysFO = 59.281^/5 - cation) and bulk 
C2/(>^ BiMnOa (Vbmo = 59AlA^/B - cation). 

The anisotropic epitaxial strain dependence of the en- 
ergies of the c-R3c and c-I4m structures is presented in 
Fig. [Ml At 0% strain, the energy difference between 
FeAFMMnFM to G-AFM in the c-R3c structure is 5 
mcV/B-cation; this differs slightly from the value for 
the relaxed structures reported in Table |Vl] due to the 
difference in lattice constants between the epitaxial con- 
strained structure and the fully relaxed structure (the 
corresponding energy difference at 0% strain for c-I4cm 
is greater because the difference in lattice constants is 
greater). At 3% strain, there is a magnetic transition 
from ferrimagnetic FeAFMMnFM to this low-lying G- 
AFM phase, while the structure remains c-R3c. This 
arises from the modification of the exchange couplings 
by the structural changes produced by the changing epi- 
taxial strain. 

With a further increase in strain, there is a transition 
from c-R3c to a c-I4cm phase at about 4.5%. Since the 
favored magnetic ordering is G-AFM in both structural 
phases, no magnetic transition accompanies the struc- 
tural transition. However, this result does illustrate the 
feasibility of a epitaxial-strain-induced structural tran- 
sition from one pattern of rotational distortions to an- 
other in this system; a coupled structural-magnetic phase 
boundary thus may be brought to light by future explo- 
ration of various epitaxial strain constraints. 



VIII. ROLE OF THE B-SITE CATION 
ORDERING IN MAGNETOSTRUCTURAL 
EFFECT 

Quasidegenerate magnetic states in the BiFcOs- 
BiMnOa nanocheckerboard are a necessary ingredient 
for the observed magnetostructural effect (cf. Fig. [T3] 
and Table IVI[) , where change in the magnetic ordering is 
achieved by a perturbation (e.g., epitaxial strain, as in 
Sec. IVIIj ). Here, wc investigate the role of the cation- 
ordering geometry in determining the spectrum of mag- 
netic states; in particular, this will show whether the 
quasidcnerature spectrum is unique to the checkerboard 
geometry. 

Magnetic energies are computed in the ideal per- 
ovskite structure of five systems: the bulk parent BiFeOa 
and BiMnOa, the BiFeOs-BiMnOs nanocheckerboard, 
the BiFcOa-BiMnOa layered (001) supcrlattice, and the 
BiFcOa-BiMnOa rocksalt structure with Fe and Mn alter- 
nating in every other unit cell ((111) superlattice) . The 
results for the checkerboard, the rocksalt, and the layered 
superlattice are presented in Table IVlIl Bulk BiFeOa and 
bulk BiMnOa exhibit behavior similar to the (001) lay- 
ered superlattice^. In these three systems, the difference 
in energy between the magnetic ground state (G-AFM in 
bulk BiFcOa, FM in bulk BiMnOa, and FeAFMMnFM 
in the (001) layered supcrlattice) and the first alterna- 
tive state is in the range 0.10 — O.lAeV/B — cation. This 
spectral gap is sufficiently large that structural changes 
cannot lower the energy of an alternative state below that 
of the original magnetic ground state. 

Indeed, for all structures considered the lowest mag- 
netic state in bulk BiFeOa and bulk BiMnOa is G- 
AFM and FM, respectively (see Figs. [HI [H])- In the 
(001) superlattice, we calculate magnetic energies for 
the G-AFM and FeAFMMnFM magnetic states in two 
structural distortions: For 1-I4cm (see Table H]), we 
find AE = -O.bOAeV/B - cation for G-AFM and 
AE = -0.553eV/B - cation for FeAFMMnFM with re- 
spect to the FeAFMMnFM magnetic state in the ideal 
perovskite cell (see Table IVII[) . For 1-R3c, wc find 
AE ^ ~0.752eV/B - cation for G-AFM and AE = 
-0.761eV/B -cation for FeAFMMnFM. For both struc- 
tural distortions considered, the lowest energy magnetic 
ordering is FeAFMMnFM. 

In contrast, all magnetic states in the nanochecker- 
board are quasidegenerate, all are lower in energy than 
the lowest-energy states in the (001) superlattice and the 
bulk. The rocksalt structure is an intermediate case: 
while the difference between the FMFM magnetic ground 
state and the first alternative state is 0.05eV/ B ~ cation, 
close to half of the spectral gap of the (001) superlat- 
tice and bulk, all the states considered fall in the same 
low-energy window as for the checkerboard. Therefore it 
is much more likely that a structurally-driven transition 
between the different magnetic states could occur in the 
checkerboard, or in the rocksalt structure, than in the 
other geometries studied here. 
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Table VII: Calculated total magnetic energies and energy differences in an ideal perovskite structure with lattice constant 
flo = 3.839A for various magnetic states in the checkerboard, rocksalt (oxygens are relaxed to accomodate their preferable 
positions), and layered superlattice of BiFeOa-BiMnOs. The checkerboard ordering shows a quasidegenerate spectrum of 
magnetic energies, whereas the rocksalt and layered superlattice show larger gaps between the ground state (GS) and the first 
alternative magnetic state. Values of Upe = Uai„ = 5eV and Jpe ~ Jaiu ~ leV are used in the first, and Up^J — AeV , 
Ulii = 5.2el/ with = f/ - J in the second column (see Sec.[Tll). 



Checkerboard E [eV/B-cation] Rocksalt superlattice E [eV/B-cation] Layered superlattice E [eV/B-cation] 
magnetic state magnetic state magnetic state 



FeAFMMnFM 

FeAFMMnFM 
FM 
C-FIM 
FeAFMMnAFM 

G-AFM 
FeFMMnAFM 



-35.04, -34.68 
AS [eV/B-cation] 
0.000, 0.000 
0.022, 0.028 
0.076, 0.113 
0.081, 0.084 
0.114, 0.152 
0.119, 0.129 



FMFM 

FMFM 
FeAFMMnFM 
FeFMMnAFM 
FM 
AFMAFM 
G-AFM 



-35.06, -34.66 
A£; [eV/B-cation] 
0.000, 0.000 
0.044, 0.070 
0.045, 0.055 
0.065, 0.047 
0.101, 0.084 
0.114, 0.079 



FeAFMMnFM 

FeAFMMnFM 

FM 
FeFMMnFM 
FeAFMMnAFM 

G-AFM 
FeFMMnAFM 



-35.11, -34.76 
AE [eV/B-cation] 
0.000, 0.000 
0.111, 0.097 
0.136, 0.143 
0.135, 0.137 
0.181, 0.219 
0.260, 0.257 



The importance of the B-site cation geometry in the 
magnetic ordering energy spectrum can be qualitatively 
understood from a simple Heisenberg model of the form 
given in Eq. (|4]), where we assume that the exchange 
couplings Ji?e, Jmji a-nd Jint arc independent of cation 
geometry, thus being transferable from one geometry to 
the other. We can approximately reproduce the mag- 
netic ordering energies in the ideal perovskite structure 
of the checkerboard with the Mn-Mn interaction Jmu be- 
ing strongly FM, the Fe-Fe interaction Jpe being AFM 
and about half the strength, and the Fe-Mn interaction 
Jint being weakly FM. Assuming the same values in the 
FeAFMMnFM GS of the (001) layered superlattice, the 
high and medium-strength bonds are all satisfied (happy 
in the language of frustrated magnetism) and the only 
unhappy bonds are weak bonds between the Mn and the 
opposite spin Fe in the adjacent layer (one bond per B 
cation). Thus this state is energetically clearly preferred 
over other ordcrings considered, which all involve a signif- 
icant fraction of unhappy high and/or medium strength 
bonds, thus opening the observed gap in the magnetic 
energy spectrum. In contrast, in the checkerboard, the 
total fraction of high and medium-strength bonds is half 
that in the layered superlattice, and the alternative states 
are low in energy as they involve tradeoffs between a 
larger number of happy weak bonds and a smaller num- 
ber of unhappy medium or high-strength bonds. Finally, 
in the rocksalt structure, all the nearest neighbor bonds 
are weak. This is consistent with the fact that all or- 
derings considered are at low energies. However, a sim- 
ple one-parameter model does not correctly account for 
the energetic order of the states in this range or the gap 
between the ground state ordering and the first alterna- 
tive state, which would require a model including next- 
nearest neighbor interactions. 

Indeed, the assumption of exact transferability used 



above is only semi-quantitatively valid. In particular, 
changes in B-site cation geometry result in relative energy 
shifts of the Fe, Mn and O states and changes in the or- 
bital wavefunctions, and thus in changes to the wavefunc- 
tion overlaps and energy denominators that contribute to 
superexchange. This leads to different values of the mag- 
netic couplings Jpe, Jmti, or Jint in the various geome- 
tries considered (cf. Table IVIip . In addition, structural 
distortions modify these magnetic couplings, as would 
be needed to explain the difference in the ordering of the 
magnetic energies in Tables IVII and IVIII However, the 
simple model does serve to give useful insight into this 
complex issue, and highlights the fact that the magnetic 
ordering spectrum is indeed very sensitive to the B-site 
cation arrangement. 



IX. DISCUSSION 

The experimental realization of the BiFcOs-BiMnOs 
nanocheckerboard would be challenging as its forma- 
tion energy is positive: the combined total energies of 
R3c G-AFM ground state of bulk BiFeOs {E[BFO] ^ 
— 35. 079eV/ B — cation) and of the R3c FM lowest energy 
structure of bulk BiMnOg {E[BMO] = -36.676eV/B - 
cation) arc lower than that of the c-R3c FeAFMMnFM 
ground state of the BiFeOa-BiMnOa nanocheckerboard 
{E[BFMO] = -71.6MeV/2B ~ cations). Though the 
BiFeOa-BiMnOa nanocheckerboard is at best metastable, 
there is indication from experiments that fabrication of 
the BiFeOa-BiMnOa nanocheckerboard with square sizes 
on the order of a unit cell would not be impossible with 
appropriate tuning of growth parameters. Growth of 
(001) BiMnOa on BiFeOa films has recently been re- 
ported. In this study post-annealing led to intermixing 
of the Fc and Mn, with a concomitant increase in ferro- 
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magnetic T^. This experiment provides support for the 
first-principles observation that magnetic ordering in this 
system is very sensitive to the B-site cation arrangement. 
With regard to other fihn orientations, (110) and (111) as 
weU as (001) BiFeOa fihns have been successfuUy grown 
on oriented SrTiOa substratea24ia. For BiMnOs, (111) 
and (001) oriented fihns have been grown with substrate 
vicinahty^^i^i. There should be no fundamental obstacle 
to analgous growth of (110) oriented films of BiMnOs. 
More generally, a combination of patterned substration, 
possible masking, layer-by-layer growth, and carefully 
tuned growth parameters could influence the deposition 
process enough to produce a checkerboard structure of 
BiPeOa-BiMnOa. 

In order to make better contact with future experi- 
ment, it is useful to consider magnetic ordering of larger- 
scale nxn BiFeOs-BiMnOs checkerboards, where the lat- 
eral dimension of the BiFeOs and BiMnOs pillars is n 
perovskite lattice constants. Within each pillar, BiFcOa 
and BiMnOa regions should be G-AFM and FM respec- 
tively, since this ordering is the most energetically favor- 
able in the parent bulk structures. This is true even in the 
extreme case of n=l discussed in Sec. |Vl The magnetic 
coupling constants Jpe, Jmu and Jmt (Eq. (O), obtained 
in Sec. IV Bl would allow the construction of Heisenberg 
models to explore the magnetic ordering of these larger- 
scale checkerboards with the ideal perovskite structure, 
where Jp,, and Jmu connecting two atoms in the xy plane 
can be taken to be equal to the coupling along the pil- 
lar, but we have not pursued this farther here. Generally 
speaking, we expect that the possibility of a structurally- 
driven magnetic transition should decrease as the lateral 
size of the BiFeOa and BiMnOa pillars increases and the 
interface effects (Fe-Mn interactions) become less impor- 
tant. 

On a technical note, the robustness of our calcu- 
lated first-principles results has been checked by using 
two different implementations of LSDA+U with different 
parametrizations to compute magnetic ordering energies 
in the ideal perovskite structure for the checkerboard, 
the rocksalt cation ordering, and the layered superlattice 
(Table VI) . The key results are the same for both imple- 
mentations: the type of ground state magnetic ordering 
for each cation arrangement, the quasidegeneracy of the 
spectrum of magnetic energies in the checkerboard, the 
gap in the energy spectrum in the layered (001) super- 
lattice, and the intermediate character of the rocksalt 



ordering. 

Finally we remark that our first principles calculations 
do not include spin-orbit coupling (SOC) that is known 
to lead to weak ferromagnetism in BiFcOgr^. Since the 
BiFeOs-BiMnOa nanocheckerboard already has a ferri- 
magnetic ground state without SOC, inclusion of SOC 
may result in a slightly changed value of the total mag- 
netization and to small canting angles of the Fe and 
Mn spins; these changes should not fundamentally affect 
the results presented here. The addition of SOC to our 
present calculations is certainly worth pursuing in future 
work. 



X. SUMMARY 

In this paper, the structure and properties of an 
atomic-scale BiFeOa-BiMnOa checkerboard were inves- 
tigated using first-principles calculations and magnetic 
modeling. This unusual heterostructurc was found to 
have properties distinct from those of its bulk parent 
constituents, or those of (001) superlattices of these two 
materials. We attribute this behavior to the magnetic 
frustration resulting from its B-sitc cation geometry; this 
leads to a quasidegenerate manifold of magnetic states 
that can be switched through small applied external per- 
turbations, resulting in an unusual magnetostructural ef- 
fect. The possibility of realizing this system in the labo- 
ratory was discussed. This study of a two-component 
nanocheckerboard should be considered as a proof-of- 
principle example, and we plan to study similar geome- 
tries on longer length scales to facilitate contact with 
future experiments. 
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